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Abstract 

We calculate the invariant-mass distribution for the cc pair produced in the inclusive T(15) decay 
based on the color-singlet mechanism of the nonrelativistic quantum chromo dynamics factorization 
approach at leading order in the bottom-quark velocity v^, in the meson rest frame. As the short- 
distance processes, we consider bh — > g*gg followed by g* — > cc and bb ^ j* ^ cc at leading 
order in the strong coupling. The invariant-mass distribution of the bb — > ccgg contribution has 
a sharp peak just above the threshold and that of the bb ^ 'j* ^ cc channel is concentrated at 
the maximally allowed kinematic end point. We predict that r[T(lS') — > cc + X]/r[T(lS) 
light hadrons] = 0.065 a^, which is smaller than a previous result by about 20%. 

PACS numbers: 12.38.-t, 13.20.Gd, 14.40.Gx 
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I. INTRODUCTION 



Recent experimental analyses being carried out by the CLEO III experiment on the 
inclusive charm production in bottomonia decays [1] have activated a series of theoretical 
studies [2-4] based on the nonrelativistic quantum chromodynamics (NRQCD) factorization 
approach [5, 6]. Predictions for the branching fractions and the charmed-hadron momentum 
distributions in inclusive Xbj {J = 0, 1, and 2) decays were presented in Ref. [2]. These 
processes are nice probes to investigate the color-octet mechanism of the NRQCD factoriza- 
tion approach which is distinguished from the color-singlet-model calculation in Ref. [7] . In 
addition to the spin-triplet P-wave (^-Pj) bottomonium decay, the CLEO Collaboration also 
analyzes the charm production in the spin-triplet S-wave {^Si) bottomonium decay. Early 
theoretical studies on the inclusive charm production in the inclusive T(15') decay began 
in the late 1970s. In 1978, Fritzsch and Streng calculated the invariant-mass distribution 
of the cc pair produced in the inclusive T(15') decay [8] by considering the QCD process 
bb — > g*gg — > ccgg within the color-singlet model, where they predicted the branching frac- 
tion Br[T(15') cc + X] to be a few percents in the limit that the charm-quark momentum 
p* in the cc rest frame can be neglected. In 1996, Cheung, Keung, and Yuan calculated 
the J/tp production rate with the same short-distance process, where they considered the 
decay of the color-singlet bb pair producing a cc pair in the color-octet spin-triplet state 
that evolves into the J/iJj [9]. One can find more studies on the bottomonium decay in 
Refs. [10-12]. 

In this paper, based on the color-singlet mechanism of the NRQCD factorization ap- 
proach, we compute the invariant-mass distribution of the cc pair in T(15') — > cc + X. This 
is an extension of a recent work on the total charm production rate and the momentum dis- 
tribution of the charm hadrons produced in the inclusive T(nS) decay [3]. As is studied in 
Ref. [3], we consider the decay of the color-singlet spin-triplet bb pair [bbi{^Si)] into g*gg fol- 
lowed by g* — > cc, which we call the QCD contribution. We also consider the short-distance 
process bbi{^Si) — > 7* — > cc, which we call the QED contribution. The perturbative calcula- 
tions of the short-distance processes are carried out at leading order in the strong coupling ag 
and the QED coupling a. For the long-distance part of the NRQCD factorization formula, we 
consider the leading contribution with respect to the bottom-quark velocity in the T(liS') 
rest frame. The relevant NRQCD matrix element for the channel is {T{lS)\Oi{^Si)\T{lS)) 
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defined in Ref. [6], where Oi{^Si) is the color-singlet spin-triplet four-quark operator for 
the annihilation decay of the T(15'). At higher orders in Vb, the color-octet short-distance 
processes bbgi^Si) g* ^ cc, bbgi^SQ) g*g ccg, and bb8{^Pj) g*g — > ccg can also 
contribute to the charm production in the inclusive T(IS') decay. The color-octet contribu- 
tions are estimated to be about 10% of the color-singlet contributions bbi{^Si) — > ccgg and 
bbi{^Si) J* ^ cc [3]. Therefore, we neglect the color-octet processes in this work. 

This paper is organized as follows: In Sec. II, we present the NRQCD factorization 
formula for the inclusive charm production in T(15') decay and compute the short-distance 
coefficients for the cc invariant-mass distribution of the process. Numerical analysis for the 
invariant-mass distribution is given in Sec. Ill, which is followed by a summary in Sec. IV. 



II. CHARM-QUARK PRODUCTION IN T(IS') DECAY 

In this section, we present the NRQCD factorization formula for the cc invariant-mass 
distribution in T(15') — > cc + X within the color-singlet mechanism at leading order in Vb. 
As the short-distance contributions, we consider the QCD process bbi{^Si) — > g*gg followed 
by g* cc and the QED process bbi('^Si) — * 7* ^ cc. We follow the formahsm to calculate 
the inclusive charm production rate in the T{nS) decay in Ref. [3]. 



A. NRQCD factorization formula 

Within the color-singlet mechanism of the NRQCD factorization formalism the differential 
rate for the inclusive decay T(15') — > cc + X at leading order in Vb is given by [3] 

dr[T(lS) - cc- + X]^ dC[^^ {0,eS,))ri,s) 

where (Oi(^5i))t(i5) = {T{lS)\Oi{^Si)\T{lS)) is the leading-order color-singlet NRQCD 
matrix element for the T(15') and rUb is the mass of the bottom quark. The short-distance 
coefficient dC['^^ is insensitive to the long-distance nature of the T(15') and calculable as a 
perturbative series of the strong coupling a^. At leading order in and a, the dimensionless 
short- distance coefficient dC[''^ is the sum of the QCD and QED contributions: 

dC}^) = dC}Q°°) + dCi'^^°). (2) 
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As is mentioned earlier, the QCD process denotes bbi{^Si) — > ccgg followed by g* — > cc and 
the QED process stands for bbi{^Si) ^ 7* ^ cc. We use the identifiers (QCD) and (QED) 
to denote those contributions, respectively, in the remainder of this paper. 

The short-distance coefficients can be determined by perturbative matching. If we re- 
place the heavy-quarkonium state |T(15')) in Eq. (1) with the perturbative bbi{^Si) state, 
\bbi{^Si)), then the long-distance matrix element {Oi{^Si))r{is) is replaced by the pertur- 
bative NRQCD matrix element {Oi(^Si))hb^(3Si), while the short-distance coefficient dC^"^^ 
remains the same: 

Both the left side of Eq. (3) and the matrix element {Oi{^Si))jji^(^3Si) calculable pertur- 
batively. By taking the ratio of the two quantities, we can determine the short-distance 
coefficient dC['^\ 



B. Short-distance coefficients 



In this section, we describe the procedure to compute the short-distance coefficient 
^^(QCD) ^j^^ QCD process. We also provide the expression for dC[^^^^ by quoting a 
previous result in Ref . [3] . In both short-distance processes, the momenta of the b and the 
b can be expressed in terms of the total momentum P and the relative momentum q of the 
bb pair: 

p ^ ^P + q, (4a) 
p=lP-q, (4b) 

where p and p satisfy the on-shell conditions p^ — p^ — ml and P ■ q — 0. In the rest frame 



of the bb pair, P = {2Eb, 0) and q = (0, q), where Eb = a/w^ + q"^. 

At leading order in ag, the full QCD amplitude for the short-distance process b{p)b{p) 
c{p')c{p')g{p2)g{p3) is given by 



= (^^^(p')r"7A.(p')4;(P2)6r.(P3) 



X 

perm 



p- P2- P3 - rrib p- p^-rrib 



u{p), (5) 
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where X^perm denotes the summation over the permutations of the three gluons attached to 
the bottom-quark hne, is a generator of the SU(3)-color in the fundamental representa- 
tion, a, b, and c are color indices for the gluons, and 62 and £3 are polarization four-vectors 
of the external gluons with momenta p2 and ^3, respectively. 

In order to extract the bbi{^Si) contribution from the bb amplitude (5), wc employ the 
covariant-projection method [13-15], which replaces the outer product u{p)v{p) of the bb 
spinors with the direct product of the color-singlet projector tti and the spin-triplet projector 
6-113, which are defined by 

TTi = -^1, (6a) 

'•^3 = ^^.i. ^ {^ + m,){]/>+2E,)^ip-m,), (6b) 
4:V2Eb{Eb + mb) 

where A^c = 3 is the number of colors, 1 is the unit matrix of the SU(3)-color, and e is 
the polarization four- vector of the bbii^'Si) state so that P ■ e = 0. The projectors (6) are 
normalized as Tr[7ri7r|] = 1 and Tr[(e • n3)(e ■ HsY] = 4poPo- 

Because we are to compute the cc invariant-mass distribution, it is convenient to write 
the amphtude (5) as the product of the vector current C for g* — > cc and the amputated 
amplitude for bbi{^Si) — > ggg*, where 



B'^ = ^Tr 



, ^, ^2—-^ /se • n3 , (7a) 

perm L P- P2- P3 - p- p3 - nib J 



Note that we do not include the color factor and the coupling in Eq. (7). At leading order 
in Vb, the amplitude for the QCD process bbi{^Si) — > ccgg becomes 



(8) 

q=0 



where i and j are color indices for the charm quark and the charm antiquark, respectively. 
In Eq. (8), we put g = to take the Ufe-lcading contributions so that Eb = mb, used 
Tr[T"T^T'^] = [d"''"^ + if"''"^)/A, and chose only the symmetric component of the color factor. 
At Vb — 0, which we take at leading order in Vb, the amplitude (8) is infrared (IR) finite in 
the soft limits of any external gluons. However, in the limit that the charm quark becomes 
massless, the amplitude (8) may acquire coUinear divergences. We will return to this point 
later in Sec. II D. 
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Squaring the amplitude (8) for bbi{^Si) — > ccgg, averaging over the spin-triplet states, 
and summing over the spins of the final states, we obtain the differential annihilation rate 
for the QCD process. The contribution of the QCD process bbi{^Si) — > ccgg to the left side 
of Eq. (3) is 

-P(QCD) _ (A^|-l)(A^c'-4) SttV 
661 - 3 2! ' 

where is the phase space element for the ccgg final state and the factors 1/3 and 1/2! 
are for the average over the initial spins and for the two identical particles (gluons) in the 
final state, respectively. The tensors B^'^ and C^'^ in Eq. (9) are defined by 

B"" = (10a) 

spins 
spins 

^ Tr[{f + m,)r{f-m,)Y] 

The polarizations for the external gluons and the bbi{^S-i) pair in B^" are summed as 

J2et{X)ef{X) = -5"^ for z = 2, 3, (11a) 

A 

^6"(A)6*'^(A) = -g-P + ^. (lib) 
A 

By substituting Eq. (9) to the left side of Eq. (3) and using the following value for the 
perturbative NRQCD matrix element 

(Oi(=^5i))65i(350 = 2N,(2E,)' = 8N,ml + 0(vl), (12) 
we determine the short-distance coefficient dC[^'^^^ as 

,^(QCD)_ (^V^^lKiV^TT^ d$4 .... 

For the QED process 661 (^5*1) 7* ^ cc, we quote the result given in Ref. [3], which was 
obtained by making use of the short-distance coefficient for the leptonic decay of T(15'): 

^^(QED) ^ |e2g2^^^2(2 + r)VT^Sil - Od^, (14) 

where cq is the fractional electric charge of the heavy quark for Q = c, b and the dimen- 
sionless variables r and ^ are defined by 

^ = ^. (15b) 



where nice is the invariant mass of the cc pair and rric is the charm-quark mass. The factor 
— r in Eq. (14) is the ratio of the phase space for the cc final state to the massless 
two-body phase space. 

C. Phase-space integral for the QCD process 

In order to compute the cc invariant-mass distribution of the QCD process, it is convenient 
to factor out the two-body phase space of the cc pair d^2 (pi p' + p') from the four-body 
phase space in Eq. (13), where pi, p', and p' are the momenta for the cc pair, the c, and 
the c, respectively. Then becomes the product of rf$2(pi — p' +p') and the three-body 
phase space (i$3(P pi +P2 +P3), where p2 and pa are the momenta for the external gluons, 
convolved with the cc invariant mass rricc as 

d<^^{P ^ p' + p' + P2 + Ps) = d^^siP ^ Pi + P2 + P3)^r^d<^>2{pi ^ p' + p'). (16) 

ZTT 

The phase space (16) can further be simplified by using the dimensionless variable ^ in 
Eq. (15b) and the scaled energy fraction Xi in the bb rest frame: 

2P-Pi 

for i — 1, 2, and 3. Rewriting the phase space (16) in terms of ^ and Xi, we get 



dxidx2dxz .P^di \pl\dVL* 



where Vt* is the solid angle of the charm quark with the three momentum p* in the cc 
rest frame. The scaled energy fraction xz can be integrated out by using the energy delta 
function. Then the cc invariant-mass distribution is obtained as a double integral of xi and 
^2, where the physical ranges for ^, xi, and X2 are given by 

r< e <1, (19a) 
2x/e< xi <l + e, (19b) 

1 / rz — ~\ 1 



-l^-x.-^xl-Ai^j < X2 <-\2-x, + sJxl-A^y (19c) 

We observe that C^" is the only factor that has the fi* dependence in Eq. (13). The 
angular integral can easily be done if we express the momenta p' and p' in terms of pi and 



7 



the relative momentum gf' as 

P' = \Pi + (20a) 
P' = \pi - q', (20b) 

where pi = {rricc, 0) and q' = (0, p*) in the cc rest frame. After integrating out the sohd 
angle fl* analytically, we get 



/ 



Eq. (21) behaves like the massive spin-1 tensor with momentum pi with pf — m^^. Substi- 
tuting Eqs. (16) and (21) into Eq. (13), setting = Ami leading order in vj,, we simphfy 
the differential short-distance coefficient for the QCD process: 

Here, the ranges of Xi and X2 integrals are given in Eq. (19) and the integrand F{^,xi,X2) 
is given by 

where the coefficients fn{xi,X2) are 

fo{xi,X2) = xl(x^ - If + xl(x2 - 1)' + xlixs - l)^ (24a) 
/i(a;i, X2) = -8 + 20{x2 + x^) - 20{x2 + x^f + 12x2X3 + 8{x2 + Xa)^ 

-14x2a;3(a;2 + x^) + xjxl, (24b) 

f2(xi, X2) = 6 - 16(x2 + X3) + 12(x2 + xsf + 2x2X3(x2 + - 6), (24c) 

f3{xi,X2) = -4 + 8(a;2 + 0:3) + 2x2a;3, (24d) 

f4{xuX2) = 2. (24e) 

Note that wc use 3:3 = 2 — — a;2 in Eqs. (23) and (24). The QCD contribution to the 
cc invariant-mass distribution in the inclusive T{1S) decay is finally obtained by replacing 
dCj"^ in Eq. (1) with dC[^^^^ in Eq. (22). 

Our result can be compared with a previous result in Ref. [8]. Neglecting the term 
proportional to in Eq. (21), we reproduce the function p(^) in Ref. [8] analytically, 

where p(^) is defined by 

^ dT[TilS)^cc + X]/dC 

T[T{1S) ^ light hadrons]' ^ 
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In Ref. [8] the authors used the order-o;^ color-singlet contribution to T(15') — > ggg for the 
r[T(15') — > light hadrons]. Imposing these approximations, we get 

where Fii^Si) is the short-distance coefficient of T(15') — > ggg at leading order in ag and 
[6, 10, 13]: 

r[T(l^) - = F,(3^0^^^2^¥^, (27a) 

Another check on the formula (22) can be done in comparison with a previous result for 
the color-octet spin-triplet cc contribution to the inclusive J/il) production in T(IS') decay: 
Eq. (22) is to be compared with Eq. (20) of Ref. [9]. The function F(^, xi, X2) in Eq. (23) is 
equivalent to Eq. (21) of Ref. [9] up to an overall factor. After considering the differences in 
the phase space and the normalization for the states, we reproduce the results in Ref. [9] at 
the leading order in a^, Vh-, and the charm-quark velocity Vc in the J/iIj rest frame (p* = 0). 



D. Massless-charm-quark limit 

As is discussed earher in this section, the short-distance coefficient dCf''^^^ (22) for the 
QCD process is free of IR and coUinear divergences as long as the charm quark is massive. 
However, in the limit where the charm quark becomes massless, rric — > 0, dCf^^^^ may 
acquire coUinear divergences. These singularities cancel only if we include the charm-quark 
contributions to the loop corrections to the gluon wave functions in the process T(15') — > ggg 
[2]. In the remainder of this section, we check if dCf^'^^^ (22) satisfies correct coUinear 
behavior in the massless-charm-quark limit. 

As the ffist step of the check, we can take the limit ^ — > on the function F{^,xi,X2) 
in Eq. (23). We find that F(0, Xi,X2) = fo{x\,X2) and this value is proportional to the 
color-singlet short-distance coefficient Fi(^S'i) in Eq. (27) for the decay T(15') ggg at 
leading order in ag and Vb [6, 10, 13]: 

Fi{'S,) = ^^'~^^^'~^^ ^lhny dx,dx2F{0,x,,X2). (28) 
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Next, we investigate the asymptotic behavior of Eq. (13) in the hmit nice 2mc and 
nic — > 0. As nice 0, approaches niccjl and gauge invariance requires pi to be 
orthogonal to B^^y. By making use of Eqs. (16), (21), and (13), we get 

B„./^*.(P-p' + p')C--^i (29) 

Note that, in this hmit, —B^^vg*^^ becomes the squared amphtude for hb\{^S\) — > ggg. Sub- 
stituting Eq. (16) into Eq. (13), and using the hmiting value (29), we find that the short- 
distance coefficient C^'^^^ is divergent logarithmically in the limit nice 2mc and nic 0. 



^(QCD) OLs 
2^ 



^ (30) 

(2mc)" 



2 m; 



where the factor 8A^c comes from the perturbative NRQCD matrix element in Eq. (12), the 
second and the third factors in the brackets are the color factor and the coupling for the 
process T(IS') ggg, respectively, 1/3! is the symmetry factor for the three gluons, and 
the factor 1/3 is for the average over the initial spin states. The quantity inside the square 
brackets in Eq. (30) is independent of nice ^-nd is finite. Simplifying the leading divergent 
term in Eq. (30), we obtain the asymptotic form of C^'^^\ Because the coUinear divergence 
is absent in the QED contribution, the colhnear divergent contribution in the QCD process 
is the same as that in the short-distance coefficient C^^\ 

Cf^-^Fi('5i)log^. (31) 

TT nic 

It is explicit in Eq. (31) that the collinear divergent contribution, which is of order a^, is 
proportional to the short-distance coefficient for T(15') ggg at order a\. 

The only order- contributions to T(15') ^ light hadrons that depend on nic except for 
the qqgg final state are the virtual charm-quark loop corrections to the gluon wave functions 
in T(15') — > ggg. The leading divergent term of the virtual correction is 

^(Cvirtual) ^ _3,n(0)Fi(35i) 

= ^Fi(^50log^, (32) 

TT /X 

where n(0) is the virtual charm-quark loop contribution to the vacuum polarization for an 
on-shell gluon and is the renormalization scale. We find that the collinear divergence 
cancels in the sum of the right sides of Eqs. (31) and (32). Therefore, + cj^'^irtuai)^ 
which is the complete nic dependent contributions to the hadronic decay of the T(IS') at 
order q;^, is free of collinear divergence. 
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III. NUMERICAL ANALYSIS 



In this section, we provide a phenomenological prediction for the cc invariant-mass distri- 
bution in inclusive T(15') decay by making use of the NRQCD factorization formula obtained 
in Sec. II. As shown in Eq. (2) , the short-distance coefficient dC^^ is the sum of the QCD and 
the QED contributions in Eqs. (22) and (14), respectively. Substituting the dCf' into the 
NRQCD factorization formula (1), wc obtain the differential rate depending on the scaled 
invariant mass ^ defined in Eq. (15b). The resultant cc invariant- mass distribution is 

—vniis) ^cc + x]^ 1 dc^'^ (g'i(^^i))T(is) .33. 

dml^ m|(^5) di ml 



where we use — rn\^^g-^ and m^^is) is the mass of the T(IS'). In our numerical analysis, we 
use the same input parameters as those used in Ref. [3], where the momentum distribution 
of the charm quark produced in the inclusive T(15') decay is studied. 

The short-distance coefficient dC^^ depends on the strong coupling ag and the ratio r 
defined in Eq. (15a). For the strong coupling, we take the running coupling Q;s(mx(is)/2) = 
0.215. As shown in Eq. (19a), the threshold of the phase space is determined by the ratio 
r = ml/m^. In order to make the end points of the phase space fit to the physical ones, we 
use mc — mo and m^ — mT(i5)/2 in evaluating the ratio r, where mo — 1-87 GeV is the 
average mass of the D° and and mr{\s) — 9.46 GeV [20] . Then the numerical value for 
the ratio becomes r = 4m|,/mY(i5) ~ 0.1563. This choice of mc and m^ for the ratio r seems 
reasonable for the open-charm production in the T(IS') decay. For the bottom-quark mass 
mfe that appears in the NRQCD factorization formulas (1) and (33), we use the one-loop 
pole mass m?, = 4.6 ±0.1 CeV. The numerical value for the long-distance NRQCD matrix 
element in Eqs. (1) and (33) is quoted from Ref. [3]: 

(0i('-Si))t(i5) = 3.07l°:2^ GeV^. (34) 

For more details of the determination of the NRQCD matrix element in Eq. (34), we refer 
the readers to Refs. [16-19]. 

The cc invariant-mass distribution from the QCD process is obtained by substituting 
Eqs. (22) and (34) into Eqs. (1) and (2), and integrating out Xi and X2 over the ranges in 
Eq. (19). We plot the QCD contribution to the invariant-mass distribution dr[T(lS') — > 
cc-l- X\/d^ in Fig. 1 as a function of the dimensionless variable ^. Note that dV\T{lS) 
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(keV) 




U ^ t3 1 

FIG. 1: Invariant-mass distribution of the cc pair produced in the inclusive T(IS') decay as a 
function of ^, which is defined in Eq. (15b). Only the contribution from the QCD process 66(^S'i) — > 
ccgg is shown. The contribution from the QED process bb{^Si) — > 7* — cc, which is proportional 
to (5(1 — is not shown. 

CC + X]/dm'^ — m^(i5^o?r[T(lS') — > cc + X]/d^. The distribution in Fig. 1 has a sharp peak 
of height 5.48 keV at ^ = 0.20, which is just above the threshold at ^ = r. The invariant-mass 
distribution of the cc from the QED process is obtained in a similar way by using Eq. (14). 
As shown in Eq. (14), the QED contribution is proportional to the delta function S{1 — ^) 
so that the contribution is concentrated at the kinematic end point = 1. Therefore, the 
contribution of the QED process is well distinguished from the QCD contribution. 

Integrating the invariant- mass distribution over the variable ^, we obtain the total pro- 
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duction rate of the cc pair in the inclusive T(IS') decay: 

r[T{lS) ^cc + X] = r(QC°) + r(Q™), (35a) 
r{QCD) ^ 1.44 ± 0.36 keV, (35b) 
p(QED) _ 2.60 ± 0.65 keV, (35c) 

where the theoretical uncertainties in Eq. (35) come from the uncertainties of m^, the 
NRQCD matrix element in Eq. (34), and uncalculatcd next-to-leading-order relativistic and 
QCD corrections, which we set to be 10 % {v^ ~ 0.1) and 21.5 % {as = 0.215) of the central 
value, respectively. The QED contribution F*^^™^ to T(IS') cc + X agrees with that in 
Ref. [3]. The QCD contribution r^^^D) ^^q^^^ ^^^^ ^^lat in Ref. [3] by about 2% because 
we omit the contribution of the process bbi{^Si) — > g*g^ followed by g* — > cc while this tiny 
contribution is included in Ref. [3]. As shown in Eq. (35), the QED contribution F^^™) 
occupies about 60% of the total charm production rate of the inclusive T{1S) decay. As is 
discussed in Ref. [3], we expect that the QCD next-to-lcading-order corrections to the QED 
process may modify the shape of the invariant-mass distribution. 

Our result for F^^^^^ can be compared with another previous result [8]. In Ref. [8], 
the ratio of the total charm production rate in the inclusive T(15') decay to F[T(15') — > 
light hadrons] is predicted to be F[T(1,5) cc + X]/r[r (IS) light hadrons] ^ 0.081 a^. 
This is greater by about 20 % than our result, 

F[T(lg) ^cc + X] 
F[T(15) ^ hght hadrons] = ^^^^ 

As we discussed in Sec. II, the authors of Ref. [8] made the approximation of neglecting 
|p*P/m^g in computing the function p(^) in comparison with Eq. (21) of this paper. The 
approximation leads to the overestimation of the height of the peak in the invariant-mass 
distribution and this is the reason for the discrepancy between the two results. To check this 
point explicitly, we carry out the same calculation with the approximations that were used 
in Ref. [8] . The calculation shows that the height of the peak increases by about 8 % from 
5.48 keV to 5.93 keV at the same horizontal position at ^ = 0.20. The difference increases as 
^ increases ranging up to 30%. As a result, the total rate without approximation is smaller 
by about 20 %. 
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IV. SUMMARY 



We have calculated the invariant-mass distribution of the cc pair produced in the inclu- 
sive decay of the T(15') based on the color-singlet mechanism of the NRQCD factorization 
formalism at leading order in the bottom-quark velocity Vb in the meson rest frame. As 
the short-distance processes, we considered the QCD process bbi(^Si) g*gg followed by 
g* — > cc at leading order in and the QED process bbi {^Si) —>■ j* ^ cc at leading order in 
a and a^. 

The QCD contribution to the cc invariant-mass distribution has a sharp peak just above 
the threshold, and that of the QED process is concentrated at the maximally allowed kine- 
matic end point. In comparison with a previous analysis on the QCD process, our prediction 
for the peak of the QCD contribution is lower than that in Ref. [8] and the total production 
rate of the cc pair in the inclusive T(IS') is smaller than that in Ref. [8] by about 20 %. The 
main reason for the discrepancy is that, in Ref. [8], the authors made an approximation of 
neglecting |PcP/^cc while we keep the full expression in Eq. (21) of this paper. 

We also investigate the coUinear divergences of the decay rate in the massless charm- 
quark limit mc — > 0. Although the decay rate of leading order in is free of both IR and 
coUinear divergences if the charm quark is massive, the rate acquires coUinear divergences 
in the limit ttIc — > 0. We have confirmed that our analytic expression for the differential 
decay rate reproduces the correct coUinear behavior in this limit so that the divergence 
exactly cancels that of the charm-quark loop corrections to the gluon wave function for the 
T(IS') — > ggg process. The sum of the two contributions are the rric dependent contribution 
to the inclusive T(15') decay into light hadrons at order a^. 
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